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ABSTRACT 



The properties of the orbit and the donor star in the high mass X-ray binary 
microquasar LS 5039 indicate that accretion processes should mainly occur via 
>■ ■ a radiatively driven wind. In such a scenario, significant X-ray variability would 

be expected due to the eccentricity of the orbit. The source has been observed at 



' X-rays by several missions, although with a poor coverage that prevents to reach 

Q ■ any conclusion about orbital variability. Therefore, we conducted RossiXTE 

^ ■ observations of the microquasar system LS 5039 covering a full orbital period 

■ of 4 days. Individual observations are well fitted with an absorbed power-law 

Q^' plus a Gaussian at 6.7 keV, to account for iron line emission that is probably 

Q ■ a diffuse background feature. In addition, we have taken into account that the 

^ • continuum is also affected by significant diffuse background contamination. Our 

• results show moderate power-law flux variations on timescales of days, as well as 
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the presence of miniflares on shorter timescales. The new orbital ephemeris of 
the system recently obtained by Casares et al. have allowed us to show, for the 
first time, that an increase of emission is seen close to the periastron passage, as 
expected in an accretion scenario. Moreover, the detected orbital variability is 
a factor of ~4 smaller than the one expected by using a simple wind accretion 
model, and we suggest that an accretion disk around the compact object could 
be responsible for this discrepancy. On the other hand, significant changes in the 
photon index are also observed clearly anti-correlated with the fiux variations. 
We interpret the overall X-ray spectral characteristics of LS 5039 in the context 
of X-ray radiation produced by inverse Compton and/or synchrotron processes 
in the jet of this microquasar. 

Subject headings: accretion, accretion disks — binaries: close — stars: individual 
(LS 5039) ~ X-rays: binaries — X-rays: individual (RX J1826.2-1450) 

1. Introduction 

LS 5039/RX J1826.2— 1450 is a high-mass X-ray binary runaway microquasar at a dis- 
tance of ~2.9 kpc from the Earth (Paredes et al. 2000; Ribo et al. 2002). Paredes et al. 
(2000) proposed that this microquasar is the counterpart of the high-energy 7-ray source 
3EG J1824— 1514 (Hartman et al. 1999), and suggested that microquasars could be the 
counterparts of some of the unidentified EGRET sources. In addition, Torres et al. (2001) 
classified this source as variable, and Bosch-Ramon & Paredes (2004) have proposed a model 
based on an inhomogeneous jet, emitting at 7-rays via inverse Compton (IC) scattering, giv- 
ing theoretical support to the possible association between the microquasar and the EGRET 
source. 

The first radio detection was reported by Marti, Paredes, & Ribo (1998) using the 
VLA, but the discovery of the radio jets was only possible when the source was observed at 
milli-arcsecond scales with the VLBA (Paredes et al. 2000). The presence of an asymmet- 
ric and persistent radio jet has been confirmed with new MERLIN and EVN observations 
(Paredes et al. 2002). The radio emission is persistent, non-thermal and variable, although 
no strong radio outbursts or periodic variability have been detected so far (Ribo et al. 1999; 
Ribo 2002). 

In the optical band, LS 5039 appears as a bright V=11.2, 06.5V((f)) star showing little 
variability on timescales of months to years (Clark ct al. 2001). Variations of ~0.4 mag 
have been reported in the infrared {H and K bands) but no obvious mechanisms for such 
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variability have been proposed (Clark et al. 2001). Using the radial velocities reported in 
McSwain et al. (2001, 2004) and newly obtained ones, Casares et al. (2005) have determined 
the most recent values for the orbital period and the eccentricity for this binary system, being 
Porb = 3.9060±0.0002 d (with to = HJD 2,451,943.09±0.10) and e - 0.35±0.04, respectively, 
with the periastron passage at phase 0.0. 

In the X-ray domain the source has been observed several times and using different 
satellites, presenting fluxes spanning one order of magnitude in the range ~ 5-50 x 10~^^ 
erg cm^^ s^^ (see Sect. 2). These fluxes imply luminosities in the range 0.5-5 x lO'^*^ erg s^^, 
assuming a distance of 2.9 kpc. Variations in the photon index have also been observed. So 
far, it has been unclear whether these variations in flux and spectral slope were related to 
long-term variations in the mass-loss rate of the companion or to accretion changes along 
the eccentric orbit. 

To determine the possible flux and photon index variations along a whole orbital pe- 
riod, we proposed X-ray observations of LS 5039 to be performed with the RXTE sateUite 
during 4 consecutive days in 2003 July. Photometric and spectroscopic optical observations 
were performed simultaneously to the X-ray ones. Regarding the photometry, Marti et al. 
(2004) have reported that no photometric variabihty above ±0.01 mag took place during six 
consecutive nights. Regarding the spectroscopy, Casares et al. (2005) have found a nearly 
constant equivalent width of the Ha absorption line of EW ~ 2.8 A, implying no significant 
changes in the mass-loss rate of the primary during our RXTE observations. 

This work has been organised as follows. In Sect. 2, we present results from previous 
X-ray observations. In Sect. 3, we show the results obtained from our RXTE data reduction 
and analysis. Further, in Sect. 4, we discuss our results in the context of both previous 
observations and the microquasar scenario and finally, in Sect. 5, we summarise this work. 

2. Previous X-ray results 

The first X-ray detection of LS 5039/RX J1826.2-1450 was made with ROSATin 1996 
October, obtaining a flux of 7.1x10"^^ erg cm^^ s~^ in the range 0.1-2.4 keV (Motch et al. 1997). 
Ribo et al. (1999) reported iiXTi? observations carried out in 1998 February 8 and 16. The 
flux detected in the energy range 3-30 kcV was ~ 50 x 10^^^ erg cm^^ s^^, and the obtained 
photon index was F =1.95, although these values are likely affected by contamination from 
the diffuse background (see discussion below). The X-ray timing analysis indicated the ab- 
sence of pulsed or periodic emission on timescales of 0.02-2000 seconds. The source spectrum 
was well represented by a power-law model plus a Gaussian component describing a strong 
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iron line at 6.7 keV. Significant emission was seen up to 30 keV, and no exponential cut-off 
at high energy was required. Unpublished observations performed by ASCA (1999 October 
4) present a photon index of about 1.6, and an unabsorbed flux in the range 0.3-10 keV of 
about 13 X 10^-*^^ erg cm~^ s^^ (Martocchia, Motch, & Negueruela 2005). iJe^^o^ylX obser- 
vations, on 2000 October 8, found a flux of ~ 5 x 10"^^ erg cm~^ s~^ (0.3-10 kcV), a photon 
index of 1.8 and a hydrogen column density of A^h = l-Olo;3 cm~^ (Reig et al. 2003a). These 
authors also concluded that an X-ray eclipse by the companion star was ruled out by these 
data. The new ephemeris obtained by Casares et al. (2005) confirm that the BeppoSAX 
observations took place during the inferior conjunction of the primary. The lack of eclipses 
provides an upper limit for the inchnation of the system of i < 69±1° (see discussion in Reig 
et al. 2003a for details). We have analysed unpublished Chandra data (2002 September 10) 
and have found a flux of 8.9 x 10^^^ erg cm~^ s~^ and F =1.14±0.12 in the range 0.3-10 keV. 
Finally, XMM observations of LS 5039 were carried out recently (2003 March 8 and 27), and 
showed that the source presented a flux of about 10 x 10~^^ erg cm~^ s~^ (0.3-10 kcV) and a 
photon index of about 1.5 (Martocchia et al. 2005). very similar fluxes and photon indexes 
are found for both observing periods, being consistent with the fact that both were taken at 
the same orbital phase (around 0.55 with the new ephemeris). 

In any but the RXTE observations in 1998, a Gaussian component in addition to the 
absorbed power-law model does not improve the fit. Therefore, it is possible that the Gaus- 
sian line of the 1998 RXTE observations could come from the Galactic ridge, being present 
because of the large field of view and effective area of the instrument. We investigate this 
issue below. For all the observations, disk features are not necessary to fit the data. 

We hst in Table 1 the dates, missions, orbital phases, photon indexes and extrapolated 
fiuxes in the energy range 3-30 keV of LS 5039 using a power-law spectrum, assuming the 
same slope all up to 30 keV, for all the prior observations but the ROS AT one (because of the 
big difference in energy range). Although it is clear that photon index and fiux changes are 
present, since the different observations are separated months/years in time, it has not been 
possible up to now to clearly establish whether the intensity of the emission changed with 
time due to variations in the mass-loss rate of the stellar companion or due to the orbital 
eccentricity. There is observational evidence supporting the first option, based on correlations 
found between the wind intensity of the companion (inferred from the equivalent width of the 
Ha absorption line) and the observed X-ray fiux (Reig et al. 2003a; McSwain et al. 2004), 
although this does not preclude the second one to occur as well. The photon index variation 
could also be a consequence of either long-term phenomena, short-term phenomena, or both. 
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3. RXTE observations 

The main goal of these observations was to perform an in depth study of the X-ray 
variabihty during a full orbital period of the microquasar LS 5039. In particular, we wanted 
to answer the following questions: Does the source flux vary with orbital phase due to a 
geometric effect related to the eccentricity of the binary system? Is the power-law index 
dependent on the orbital phase? Does the iron line come from the source? This section tries 
to answer these three queries, although a deeper and wider discussion is developed in Sect. 4. 

3.1. Data reduction and analysis 

RXTE observed the source on 17 runs of different duration covering a total of 4 days 
between 2003 July 4 and 8. Due to the relatively faint emission of LS 5039 during the 
observations, the source spectrum is background dominated beyond 30 keV, so HEXTE 
data has not been considered and we have focused on the analysis of PCA data. For the 
four longer runs, only PCU # 2^ data was available, whereas two and even three detectors 
were on during significant fractions of time of the other runs (PCUs #1, 3 and 4). All 
the observations performed with PCA have been reduced and analysed. A summary of the 
observations is shown in Table 2. 

We have reduced the data using the HEAsoft FTOOLS package version 5.3.1. The 
spectral analysis has been performed with Xspec 11.3.1. To carry out the data reduction, 
we have created a filter file to avoid data taken with an elevation angle with respect to the 
horizon smaller than 10 degrees, electron contaminated data, data taken during the South 
Atlantic Anomaly passage and data with an offset of the pointings bigger than 0.02 degrees. 
We have accounted also for the detectors that were on during the observation. 

For the 17 runs, we have tried to fit the spectra with an absorbed power-law plus 
a Gaussian emission line and a disk black body components. For all of them either the 
difference between models with and without a black body component is insignificant or the 
black body flux is consistent with being zero. Therefore, to estimate the photon index and 
the fluxes of both the continuum and the line, we have fltted the data only with an absorbed 
power-law plus a Gaussian line around 6.7 keV with fixed zero width (due to the low spectral 
resolution of RXTE). Regarding absorption, we note that inhomegeneities in the velocity 



^Concerning PCU #0, despite efforts to recalibrate its response function following the loss of one of its 
gas layers, its spectra continue to deviate significantly at low energy compared to other PCUs. Therefore, 
although PCU # was also on during the runs, we have not included it in our spectral analysis. 
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and/or spatial distribution of the stellar wind could produce a variable amount of hydrogen 
column density, as has been observed in other wind- fed accreting systems like GX 301—2 
(White & Swank 1984). Unfortunately, it is not possible to extract such kind of information 
from our RXTE data of LS 5039 due to the faint emission of the source, the relatively low 
value of A^H and the low-energy threshold of our data at 3 keV. Therefore, we have assumed 
a constant hydrogen column density and fixed it to a value of 8.7 x 10^^ cm~^ (resulting 
from the average between the one found by Reig et al. (2003a) using the BeppoSAX data 
and the one derived from the formula of Predehl & Schmitt (1995) apphed to the quoted 
value for Ay by McSwain et al. (2004)). The spectrum and the fit to the data taken at 
phase 0.78-0.84, during a long-lasting maximum of emission, are presented in Fig. 1. The 
aspect of the spectra and the fits obtained from the other observations are quite similar to 
the one shown here. The date, time and orbital phase of the 17 runs, with their power-law 
indexes and unabsorbcd fiuxes for the power-law component in the energy range 3-30 keV 
are presented in Table 2. 

Due to the low source flux, the upper limits on the source variability within a given 
run are not very constraining. The derived limits depend on a number of factors, including 
the fitting functions used (power-laws, Lorentzians, etc.) and the background subtraction 
(total background is about 30 cts s^^ PCU~^ in the 3-30 keV energy range). Fast Fourier 
transforms of the high time resolution data from all active PCUs were performed to create 
power density spectra (PDS) in the 0.01-1000 Hz range. Given the faint nature of the 
source, all available energy channels were combined when making the PDS, and the obtained 
fractional variabihty amphtude is generally less than 25% (Scr). The PDS are dominated by 
Poisson noise, and generally appear to have no significant noise components which can be 
attributed to the source. 



3.2. Results 

The unabsorbed power- law flux and the photon index as a function of time are presented 
in Fig. 2. We flnd evidence of day to day variations in both parameters. We show the same 
data folded with the orbital period in Fig. 3. It can be seen that the flux appears to be 
moderately variable along the orbit, presenting a variation of a 50% (see Table 2) with a long- 
lasting maximum at phase 0.8 and a short-duration maximum around phase 0.16. This short- 
duration maximum, of about half an hour, is more than three sigmas over the level of emission 
indicated by the adjacent points, and it is not unique since other short-duration peaks are 
present even within the longer runs. The plotted flux in Fig. 3 has been corrected for the 
internal background. However, this background subtracted flux still contains a contribution 
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from the diffuse background that is difficult to disentagle from the source spectrum and is 
estimated to be about 2 counts per second per PCU^. Applying an absorbed Raymond plus 
power-law model (see Valinia & Marshall 1998) to this count rate in the 3-30 keV energy 
interval, we obtain a diffuse background flux of about 3x10^^^ erg cm^^ s"^. Thus, the ratio 
between the long-lasting maximum and the minimum fluxes after subtracting the diffuse 
background is ~2.5. 

In the lower panel of Fig. 3, we show the evolution of the photon index with the orbital 
phase. As can be seen, the photon index varies along the orbital period, presenting a 
hardening close to pcriastron. Actually, comparing in Fig. 3 the unabsorbcd power-law 
flux and the photon index evolutions at different phases of the orbit, an anti-correlation 
between both appears. This is better illustrated in Fig. 4, where we plot the photon index 
as a function of the unabsorbed power-law flux. It seems that, when the flux is higher, 
the photon index is harder, and the contrary happens when the flux is lower (this anti- 
correlation holds even within the longer observations when split and analysed by parts) . We 
stress here that changes in flux and photon index cannot be due to changes in the hydrogen 
column density for this weakly absorbed source, since we are working in the energy range 
3-30 keV, nor to diffuse background variations (internal background ones are already taken 
into account) since, after applying the diffuse background model, the anti-correlation is still 
there. We note that the photon indexes calculated using this model are harder, in the 
range ~1.3-1.6, more in accordance with those obtained in previous X-ray observations of 
the source, excluding the 1998 RXTE ones. Nevertheless, the uncertainty in the diffuse 
background model determination at the LS 5039 position requires to be cautious about 
these results. Moreover, the flux variability does not seem to depend only on the photon 
index variability, since when flxing it to an intermediate value of 1.9 the flux variation is 
still signiflcant. Another test has been done by fltting the unabsorbed flux data given in 
Table 2 with a constant value. The obtained reduced greater than 3 indicates that the 
fit is not good. Moreover, the diffuse background can not introduce variability, since its 
value, estimated by the RXTE team, remains constant within about a 10% during periods 
of time longer than the orbital period. Finally, no significant optical variations have been 
observed in simultaneous photometric observations (Marti et al. 2004), which indicates that 
the companion star showed a stable behavior from the photometric point of view during our 
RXTE observations, reinforcing the idea that the measured X-ray fiux variability is due to 
changes of the accretion rate in an eccentric orbit. 

Although the data do not cover continuously the whole orbit, there appear to be two 



^Although this value has been found by the RXTE team after observations performed during the 2004 
fall, we have assumed that during our observations, performed in 2003 July, it was roughly the same. 
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types of variation timescales in the lightcurves. One apparently linked to the orbital period 
of ~3.9 days, the other of shorter duration, like miniflares, of about 1 hour. No significant 
changes are detected through simultaneous (although not strictly) spectroscopic observations 
of the Ha line along the orbit on day timescales (Casares et al. 2005). Thus, stellar wind 
changes can be ruled out to explain the first type of variability. In contrast, fast fluctua- 
tions in the stellar wind might be related to the shorter timescale variability as significant 
variations on timescales of minutes have been observed in radial velocity measurements of 
LS 5039 (see figure 2 of Casares et al. (2005), around HJD 2,452,827.7). Unfortunately, 
there is no spectroscopic data right before these minifiares observed in X-rays. 

Our fits to the spectra are well consistent with the presence of a line at 6.7±0.1 kcV in 
the data, as in the previous iiXTE observations (Ribo et al. 1999). However, as mentioned 
above, no other X-ray mission but RXTE (which does not have imaging capabilities) reported 
the presence of the iron line, so its origin could be in the Galactic plane regions behind the 
source, and due to the so-called Galactic ridge emission. One possible test for determining 
whether or not the line is related to LS 5039 would be the presence or not of line-fiux 
variations correlated with continuum fiux changes. In fact, the different values of the line 
fiux obtained on the different runs are always consistent, at 90%confidence, with their average 
value of ~ 1.65 x 10^"^ photon cm~^ s~^ per PCU beam of 0.9 square degree. For comparison, 
the value of the iron-line flux at the LS 5039 position after measurement with the Ginga 
Large Area Proportional Counters (LAC) is ~ 4.0 x 10"^ photon cm~^ s"^ per LAC beam 
of l°x2°, with the large beam oriented perpendicular to the Galactic plane (Yamauchi & 
Koyama 1993). Therefore, this implies a flux around 1.8 x 10~^ photon cm~^ s~^ per 0.9 
square degree, compatible with our individual measurements and slightly above the average 
value we have obtained. Moreover, we have performed the analysis of the slew data (when 
RXTE was pointing off the source) along the four days of observation and found marginal 
evidence of the presence of a line (3cr) with the flux similar to that obtained from the on- 
source data. Although slew data have been difficult to model and we had to restrict the 
energy range to 3.0-10 keV, we have fitted first the data to a Bremsstrahlung obtaining a 
reduced of 2, and afterwards to a Bremsstrahlung plus a Gaussian line, improving the 
reduced value down to 1. The former facts seem to rule out the possibility that the 
line comes from the source, although a deep observation by an instrument with imaging 
capabilities, if possible around periastron passage when the fiux is higher, is still required to 
definitively solve this question. 
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4. Discussion 

4.1. Orbital X-ray variability 

The RXTE observations presented here show evidence of flux variability occurring on 
timescales of days which is anti-correlated to the photon index of the spectrum. Since this 
happens while no changes are detected in the companion star and the maximum X-ray flux 
occurs near pcriastron, we propose that this variability is related to the orbital motion of the 
compact object accreting from the companion star wind along an eccentric orbit. Actually 
the emission peak seems to take place at phase 0.8, not right at periastron passage. Leahy 
(2002) models a similar X-ray lightcurve observed in GX 301—2 through a non spherically 
symmetric stellar wind accreted by the compact object following Bondi-Hoyle accretion. 

We can consider that the wind of the companion star follows a /5-law and use a simple 
spherically symmetric Bondi-Hoyle accretion model (see Reig et al. 2003a for details) to 
study the ratio of maximum to minimum flux, to be compared to the observed one. For 
this purpose we have assumed the following parameters: Mopt = 40 M©, i?opt = 10 Rq, 
Mopt = 1 X 10~^ Mq yr~^ (average between the value given in Kudritzki & Puis 2000 from the 
wind-momentum relationship and the value inferred from optical spectroscopy by McSwain 
et al. 2004), Voo = 2440 km s"^ (from observations by McSwain et al. 2004), (3 = 0.8, and 
flnally a neutron star of Mx = 1.4 Mq with a radius of Rx = 15 km in an eccentric orbit 
with Porb = 3.9060 days and e = 0.35. (Casares et al. 2005). These parameters provide an 
expected variability with a ratio of ~10 between the maximum and minimum flux. Even 
when considering the uncertainties in all the involved parameters, ratios between 8-17 are 
found. These ratios are much lower than the variability ratio ~40 considering the previous 
eccentricity value of about 0.5 (McSwain et al. 2004), but still too large to explain the ratio 
of ~2.5 observed by RXTE. 

A possible solution that would easily reconcile the theoretical and the observed flux 
ratio is the presence of an accretion disk around the compact object. This would smoothen 
the accretion rate of the compact object, implying lower and longer emission bumps around 
periastron passage in the hghtcurves. Therefore, the difference between the expected and 
the observed flux ratio could be the flrst hint, in the X-ray domain and independent of the 
presence of relativistic jets, of the existence of an accretion disk in LS 5039, although it is not 
clear what kind of disk can be present at such a low ratio of X-ray to Eddington luminosity. 
It should be noted that for the short-timescale variations or miniflares, although they present 
the same flux/photon index anti-correlation as the rest of the observations suggesting to be 
intrinsic to the source, the disk smoothing might not operate in the same manner as for long- 
lasting peaks, since it is the wind itself that changes rapidly. This could affect the accretion 
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process in a more violent manner (e.g. inducing unsteady accretion and subsequent sudden 
X-ray flux increase) than a slow changing flow, as it would be the case for an approximately 
constant mass-loss rate crossed by the accreting compact object. 

BeppoSAX observations show a decrease in flux of about a 50% between the periastron 
passage at phase 0.0 and the end of the observations at phase 0.21 (flgure 1 of Reig et al. 
2003a). Interestingly, this is not too far from the flux reduction in this phase interval 
predicted by the accretion model discussed above, and we remind the reader that BeppoSAX 
had imaging capabilities, preventing the influence of background sources in the data. We note 
that the position of the maximum of emission in the BeppoSAX observations is consistent to 
be before periastron, like in the observations presented here. Therefore, these observations 
are complementary to the RXTE ones reported here and give clear support to the orbital 
variability as due to the eccentric orbit. 

Concerning the long-term correlations in the X-ray emission and the stellar mass- 
loss rate, the new values for the equivalent width of the Ha absorption line of ~2.8 A 
(Casares et al. 2005) are similar to those obtained during the epochs when ASCA and 
Chandra observed the source. Moreover, taking into account the diffuse background emis- 
sion, the average X-ray flux truly coming from the source during our observations is about 
2xl0~^^ erg cm~^ s^^, similar also to the fluxes observed during ASCA and Chandra ob- 
servations. Therefore, the apparent correlation between changes in the mass-loss rate of the 
stellar companion and the X-ray emission, suggested by Reig ct al. (2003a) and supported by 
McSwain et al. (2004), seems to be conflrmed by these new spectoscopic and X-ray results. 

4.2. Proposed origin for the X-ray emission 

The signature of the existence of an accretion disk around the compact object is not 
present in any of the obtained X-ray spectra of the source and no cut-off is present in data 
up to 30 kcV. Actually, the possible detection of LS 5039 by instruments working in the 
energy range beyond the RXTE one (BATSE, Harmon et al. (2004); COMPTEL, Strong 
et al. 2001; EGRET, Paredes et al. 2000) seems to preclude the existence of any cut-off 
at energies beyond the RXTE range, typical of the thermal Comptonization spectra for 
emission produced in a corona. Therefore, we propose that the X-ray emission is due to 
inverse Compton (IC) and/or synchrotron processes within the relativistic jets of LS 5039. 
In our scenario, an increase in the accretion rate would imply a higher electron acceleration 
efficiency in the jets, filling with fresh particles the higher energy range of the electron 
energy distribution. This would harden the X-ray spectrum, which is actually observed 
by RXTE through the flux/photon index anti-correlation. The smooth evolution of the 
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flux could be explained mainly, as it is mentioned above, by the dynamical timescales of a 
faint underlying accretion disk, which would smoothen the emission variation around the 
periastron passage. It has been considered by several authors that relativistic electrons in 
the jet of a microquasar could produce significant X-ray emission by IC upscattering of the 
stellar photons or their own synchrotron photons (see Kaufman Bernado, Romero, & Mirabel 
2002; Georganopoulos, Aharonian, & Kirk 2002; Reig, Kylafis, & Giannios 2003b; Bosch- 
Ramon, Romero, & Paredes 2005). Optically thin synchrotron radiation from a jet could 
also dominate at X-rays, depending on the matter density and magnetic field conditions 
within the jet (Markoff, Falcke, & Fender 2001). Jet dominated states could be common and 
the disk radiation would appear covered by the non-thermal jet emission (Fender, Gallo, & 
Jonker 2003). These jet dominated states are associated with the low-hard state of black 
hole X-ray binaries. The nature of the compact object in LS 5039 is still unknown (although 
in could be a black-hole, see Casares et al. 2005) and its relatively high radio flux density 
and optically thin radio spectrum do not agree with the observed values found in black 
holes (Gallo, Fender, & Pooley 2003). However, the observed X-ray spectrum is always very 
similar to the ones of black holes in the low-hard state: a power-law with a photon index 
around 1.5 (see Tables 1 and 2). 



5. Summciry 

We have reported recent RXTE observations of LS 5039 covering a full orbital period of 
4 days. The obtained results show moderate day to day variability that is interpreted as due 
to changes in the accretion rate because of the orbital motion of the compact object along an 
eccentric orbit, with an increase of emission right before periastron passage. In addition, the 
observed short-timescale flux variations (miniflares) are likely related to observed fast varia- 
tions in the stellar wind properties. These observational results give support to the accretion 
scenario, making other possibilities like the young non-accreting pulsar scenario (Martocchia 
et al. 2005) to be more unlikely. In parallel with the power-law flux variability, the flux of the 
iron line, detected at ~6.7 keV like in previous RXTE observations, has remained constant 
within the errors and similar to the expected background line flux. Furthermore, slew data 
of these RXTE observations present marginal evidences of the iron fine as well, and the fine 
has not been detected by X-ray satellites with imaging capabilities. All these facts make us 
suggest it is a background feature. No disk signatures are found in the spectra, although this 
is not rare for other X-ray binaries in the low-hard state. In fact, the presence of an accretion 
disk could explain the apparent smoothcning of the flux changes along the orbit observed by 
RXTE with respect to what is predicted by using a Bondi-Hoyle accretion model. Finally, 
there appears to be an anti-correlation between the flux and the X-ray photon index along 
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the orbit, reaching a maximum of emission and the hardest spectrum around phase 0.8. 
We suggest that a scenario where the X-rays come from IC interactions and/or synchrotron 
emission within a relativistic jet could explain this anti-correlation, as well as the absence of 
cut-off at hard X-rays. Although further theoretical work must be done, new observations to 
constrain the X-ray flux variability between periastron and apastron passages are necessary 
to better understand the orbital variability of LS 5039. 
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Table 1. Summary of previous X-ray observations. 



Extrapolated 3-30 keV Flux 



Date Mission Phase Photon Index (xlO erg cm ^ s ^) 



1998 Feb 08 RXTE 0.803-0.837 1.95±0.02 5.6 

1998 Feb 08 RXTE 0.006-0.057 1.95±0.02 4.6 

1998 Feb 16 RXTE 0.039-0.059 1.95±0.02 4.6 

1999 Oct 04 ASCA 0.382-0.471 ~1.6 2.7 

2000 Oct 08 BeppoSAX 0.969-0.205 1.8±0.2 0.8 

2002 Sep 10 Chandra 0.690-0.720 1.14±0.20 2.3 

2003 Mar 08 XMM 0.526-0.557 l-56tom 2.0 
2003 Mar 27 XMM 0.533 0.561 1.19l" m 2.0 



Note. — 



Errors are within the 90% of confidence. 
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Table 2. Summary of RXTE observations. 



Middle Time 


Duration 






Unabsorbed 3-30 keV Flux 


(MJD) 


(day) 


Phase 


Photon Index 


(xlO- 


'^^ erg cm ^ s ^) 


52824 61 


0.24 


782-0 842 


1 84+^-^i 

-'-•°^-0.04 








n no 


n 888 n 8Q/1 
U.ooo— U.oy^t 


^•U4_o.o9 




5.09l°:l 


52824.99 


0.01 


0.907-0.910 


2.031°:}} 




c- 1 4+0.30 


52825.06 


0.01 


0.924-0.926 






/I 09+0.30 

^•'J^-0.41 


52825.18 


0.05 


0.950-0.963 






5.48i:°:i 


52825.60 


0.24 


0.034-0.095 


2.09i°:°l 




4.421°:}° 


52825.98 


0.01 


0.159-0.162 


-'-•"-'^-0.08 




6.951°:^^ 


52826.04 


0.02 


0.176-0.181 


1 qa+0.12 




4.541°:^^ 


52826.17 


0.05 


0.202-0.217 


1 q8+°-^2 

-■-•^0-0.12 




rr 1 4+0.33 
'-'•-'-^-0.38 


52826.69 


0.17 


0.321-0.364 






p; 1 9+0.14 
^•^^-0.23 


52826.96 


0.02 


0.411-0.415 


2.081°:}J 




4.82l°:H 


52827.03 


0.01 


0.429-0.431 


9 iq+0.11 
^•-'-'-'-0.10 




4.581°:^^ 


52827.12 


0.08 


0.445-0.465 


9 1«+0.12 




4.511°:^° 


52827.57 


0.24 


0.539-0.600 


9 Q4+0.02 




4.981°:}^ 


52827.95 


0.01 


0.663-0.667 


9 ni+0.11 
^•'-'^-0.10 




5.02l°:^« 


52828.01 


0.02 


0.680-0.685 


9 nfi+0 09 




5.09l^:| 


52828.18 


0.02 


0.724-0.728 


i-yu_o.o8 




5.80_o:24 



Note. — 



Errors are within the 90% of confidence. 
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I I 1 I I I 1 I I I 

5 10 20 

channel energy (keV) 

Fig. 1. — X-ray spectrum of LS 5039/RX J1826.2— 1450 obtained with the data taken be- 
tween phases 0.78 and 0.84. The sohd hne represents the fit to the data. The lower panel 
shows the ratio between observed and fitted fluxes. 
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Fig. 2. — Unabsorbed power-law flux in the energy range 3-30 keV (upper panel) and 
photon index (lower panel) along the observational period, without subtracting the diffuse 
background model. The vertical error bars show the 90% of confidence, the horizontal error 
bars show the duration of the observation in time. The estimated source zero level is indicated 
by a dashed line in the upper panel. 
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Fig. 3. — Same as Fig. 2 but as a function of the orbital phase, computed by using P^rb = 
3.9060 ± 0.0002 d and to = HJD 2,461,943.09 ± 0.10. The periastron passage is at phase 0.0, 
and two orbital periods are shown for a better display. 
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Fig. 4. — Photon index as a function of unabsorbed power-law flux in the energy range 
3-30 keV, without subtracting the diffuse background model. The errors are within the 90% 
of confidence. 



